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• Premise of the study: Microsatellite loci were developed for the haploid lichenized fungal species Usnea subfloridana to study 
its population subdivision and the species' response to forest disturbance, fragmentation, and environmental pollution. 

• Methods and Results: We developed 14 polymorphic microsatellite markers using 454 pyro sequencing data of U. subfloridana. 
The number of alleles per locus ranged from three to 15, and Nei' s unbiased gene diversity averaged over nine markers without 
null alleles ranged from 0.64 to 0.67. Evaluation of the cross-species amplification in U. glabrescens and U. wasmuthii indi- 
cates that these markers are also informative in other Usnea species. 

• Conclusions: These markers will allow us to investigate the effects of forest management and environmental pollution on ge- 
netic population structure of U. subfloridana and closely related species. Moreover, they will help facilitate phylogeographic 
studies of U. subfloridana across the species' distribution area in Europe. 

Key words: Ascomycetes; lichen-forming fungi; microsatellites; population subdivision; Usnea florida; Usnea 
subfloridana. 



Conspicuous epiphytic lichens like beard lichens (Usnea Dill, 
ex Adans.) are valuable indicators of forest ecosystems, hence 
contributing to monitoring the conservation value of forest land- 
scapes (Will-Wolf et al., 2002). Usnea subfloridana Stirton is 
a widely distributed species occurring in Europe (Torra and 
Randlane, 2007; Randlane et al., 2009), appearing from the north- 
ern boreal to temperate regions (Halonen et al., 1998). The sexually 
reproducing U. florida (L.) Weber ex F. H. Wigg. and U. subflori- 
dana, which has a predominantly asexual reproduction with sym- 
biotic propagules, represent a typical species pair, as they share 
many morphological characters but differ by the characters as- 
sociated with their dissimilar dispersal strategies (Articus et al., 
2002; Randlane et al., 2009). Here, we develop 14 microsatellite 
markers to study the impact of land use and habitat fragmentation 
on the lichen's dispersal and population subdivision. 

METHODS AND RESULTS 

We selected three U. subfloridana specimens sampled from Norway, Fin- 
land, and Lithuania and two U. florida specimens sampled from the United 
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Kingdom (Appendix 1). The central axis, which is of pure fungal origin (hap- 
loid), was manually separated and used for DNA extraction with the DNeasy 
Plant Mini Kit (QIAGEN, Hilden, Germany), according to the manufacturer's 
protocol. Library preparation and whole genome 454 pyrosequencing of 
pooled DNA samples was performed by Microsynth (Balgach, Switzerland) 
using a Roche GS FLX sequencer to generate enough random sequences and 
to isolate a sufficient number of microsatellite loci. Shotgun libraries were 
prepared using the GS FLX Titanium Rapid Library Preparation Kit (Roche 
Diagnostics, Basel, Switzerland), and Microsynth provided barcode adapters. 
We obtained 85,718 reads with an average length of 391 bases and a total of 
27,344,042 bases out of a 2/16th run. We screened for all sequence motifs of 
di-, tri-, tetra-, and pentanucleotide microsatellites in the unassembled reads 
using MSATCOMMANDER 1.0.2 alpha (Faircloth, 2008). Microsatellites 
with motifs repeated at least eight times (for di- and trinucleotides) or six 
times (for all others) were selected. Primer pairs were developed with Primer3 
(Rozen and Skaletsky, 2000), implemented in the software MSATCOM- 
MANDER 1.0.2 alpha using the default parameters except for the following: 
optimal primer size 20 bp, product size 150-450 bp, melting temperature (r m ) 
58-65°C. We found 132 primer pairs that fulfilled the specified primer param- 
eters, but 68 pairs were later discarded either because they were duplicates, 
which were detected after alignment using CLC DNA Workbench 5 (CLC 
bio, Aarhus, Denmark), or because they contained mononucleotide repeats in 
the flanking region. 

Additionally, we set up axenic algal cell cultures of the photobiont of U. 
subfloridana to assess the symbiont specificity of the newly designed markers. 
The culture was established under sterile conditions on 1/4 of strength of origi- 
nal Trebouxia Organic Nutrient Medium-I according to Ahmadjian (1967). 
Algal cells were taken from the algal layer of U. subfloridana thalli and inocu- 
lated on the medium. The cultivation took place under diurnal light (12 h) and 
darkness (12 h) for four months before the algal culture was harvested and de- 
posited at the Swiss Federal Research Institute WSL (cultures TTA1 and 
TTA2) at -80°C. Algal cells were disrupted and DNA was extracted with MO 
BIO PowerPlant DNA Isolation Kit (MO BIO Laboratories, Carlsbad, Califor- 
nia, USA) according to the manufacturer's protocol. The three loci that 
produced positive PCR reactions were excluded from further analyses because 
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they were considered alga-specific rather than fungus-specific. For PCR ampli- 
fication, forward primers were labeled with an Ml 3 tag (5'-TGTAAAACGAC- 
GGCCAGT-3') (Schuelke, 2000). PCR reactions were performed in a total volume 
of 10 jiL containing 1 jiL of -1-5 ng genomic DNA, 0.5 jlxL of 5 pM forward 
and reverse primers, and 2x Type-it Multiplex PCR Master Mix (QIAGEN). All 
PCRs were performed on Veriti Thermal Cyclers (Life Technologies, Carlsbad, 
California, USA). The PCR reactions were assessed using a temperature gradi- 
ent with one-degree increments from 56-61 °C, and under the following condi- 
tions: denaturation for 5 min at 95°C, followed by 33 cycles of 30 s at 95°C, 90 s 
at 56-6 1°C, and 30 s at 72°C; then for the M13-tag binding an additional 
eight cycles of 30 s at 95°C, 90 s at 53°C, and 30 s at 72°C, with a final exten- 
sion of 30 min at 60°C were run. 

Only primers failing to amplify with DNA extracted from the axenic algal 
culture were considered of fungal origin. These 61 primers were tested for 
variability under the same conditions as above and using the total DNA of 
eight specimens of U. subfloridana collected from EVO population from 
southern Estonia (Appendix 1), resulting in 14 loci with satisfactory amplifi- 
cation. Cross-species amplification of two closely related species (Saag et al., 
2011) was tested with 12 specimens of U. glabrescens (Vainio) Vainio and 14 
specimens of U. wasmuthii Rasanen, which were collected from the same 
site. Approximately 50 mg dry weight of each lichen thallus was lyophilized 
overnight and ground in a Retsch MM2000 mixer mill (Diisseldorf, Ger- 
many) for 3 min at 30 Hz, and total DNA was extracted with the same proce- 
dure as the algal cells. 

To characterize the 14 polymorphic U. subfloridana loci (Table 1), we 
analyzed PCR products of 174 individuals from three populations (Table 2, 
Appendix 1). Fluorescent forward primers were used for the PCR protocol 
and the reaction was adjusted to: 5 min at 95°C, followed by 28 cycles of 
30 s at 95°C, 90 s at 57°C, and 30 s at 72°C, with a final extension of 
80 min at 60°C. PCR products were multiplexed (Table 1) and run on a 
3130x1 DNA Analyzer with GeneScan 500 LIZ Size Standard (G5 dye set) 
for fragment analysis (both by Life Technologies). Alleles were determined 



using GeneMapper version 3.7 (Life Technologies). To characterize the 
variability of each microsatellite locus, we counted the number of alleles 
and calculated Nei's unbiased gene diversity using Arlequin version 3.11 
(Excoffier et al., 2005). 

Sequences of the 14 polymorphic loci were deposited in GenBank as they 
appear in the original pyrosequencing sample (Table 1). Five loci (Usl0-Usl4) 
had more than 10% of null alleles, possibly because of mutations in the 
primer regions, and were therefore omitted from the population analyses. 
The nine microsatellite loci that had no null alleles (Usl-Us9) produced 
three to 15 alleles per locus with a mean of 8.78. Nei's unbiased gene diver- 
sity, averaged over nine markers, ranged from 0.64 to 0.67 (Table 2). After PCR 
optimization for the annealing temperature (Table 1), all nine primers success- 
fully amplified and were polymorphic in U. subfloridana, U. wasmuthii, and U. 
glabrescens, except marker Us07, which showed no polymorphy in U. glabres- 
cens. As is often the case in populations of highly clonal organisms such as li- 
chens (Walser et al., 2004; Dal Grande et al., 2012), significant linkage 
disequilibrium was found using Arlequin version 3.11 in 41 U. subfloridana 
distinct multilocus genotypes for two pairs of markers (i.e., Us02-Us06 
and Us05-Us08). 



CONCLUSIONS 

The manual separation of the purely fungal central axis of the 
genus Usnea did not provide pure fungal DNA as expected. 
This preparation led to mixed DNA of the two fungal and algal 
symbionts and thus symbiont-specificity of genetic markers has 
to be tested in lichens (Devkota et al., 2014) even if lichens 
contain purely fungal plectenchyma. The newly developed, 
highly variable fungus-specific markers are currently being 



Table 1. Characteristics of 14 microsatellite loci developed for the lichen fungus Usnea subfloridana: 











Fluorescent 




Allele size 


Primer cone. 




GenBank 


Locus 




Primer sequences (5'-3') 


Repeat motif 


dye 


r a (°o 


range (bp) 


(MM) 


Multiplex 15 


accession no. 


UsOl 


F: 
R: 


CCCTCCAATCCCATATATAC 
GATCCATCGTCTCCATGATA 


(CA) 10 


PET 


59 


181-303 


0.02 


1 


KC912190 


Us02 


F: 
R: 


GTTTGCGAGGGATACTTTAC 
GTAGAAAGACAAGGACGCAC 


(TA) 10 


NED 


57 


83-89 


0.06 


2 


KC912191 


Us03 


F: 
R: 


CCATCTCAAACCCGATAGTC 
CCTCTTATCCTTTGGCAAAT 


(CTT) n 


PET 


57 


192-218 


0.12 


2 


KC912192 


Us04 


F: 
R: 


ATTCCTCCGCAGTCAAGTAG 
CTTTACTTTGCGAGCTTCAG 


(CAT) 9 


VIC 


57 


181-196 


0.06 


2 


KC912193 


Us05 


F: 
R: 


TCTGGGATGGATGGCATT 
GCGGGAAAGCCACAGATA 


(CA) 9 


FAM 


57 


161-175 


0.20 


2 


KC912194 


Us06 


F: 
R: 


CGATGGCGTAATGATCAG 
CAAACATTCGCGTCAGTTAC 


(GAT) 9 


FAM 


57 


106-127 


0.12 


2 


KC912195 


Us07 


F: 
R: 


TCCTGAGTTTGAGCCTCTTG 
GTCCAACGCAGCTACTCTCT 


(GTC) 8 


VIC 


57 


123-135 


0.02 


3 


KC912196 


Us08 


F: 
R: 


CCTTCACCTTCCCTATACACC 
GAATCCCAGTCAATGCACTC 


(AAC) 10 


NED 


57 


192-217 


0.03 


3 


KC912197 


Us09 


F: 
R: 


GGAAGGTATGAGGTGAGGTG 
AGAATTTGCAGCACGGGTA 


(GGAT) 6 


FAM 


57 


341-372 


0.30 


3 


KC912198 


UslO 


F: 
R: 


ACGAGGGAGGGTAGTTAGG 
GTATGATTGGTTGTGATGCG 


(GA) 9 


PET 


59 


220-318 


1.0 


1 


KJ446958 


Usll 


F: 
R: 


AGCATCCCACAACACACTTC 
CACTCTTTATCCTGCTCGC 


(GCA) 9 


NED 


59 


242-245 


1.0 


1 


KJ446959 


Usl2 


F: 
R: 


CCCTTCTTCCTACTCGGTATC 
CCTACTTCTACCATGCGTCC 


(TTC) 8 


PET 


57 


281-287 


0.2 


3 


KJ446960 


Usl3 


F: 
R: 


CCATTTACGGCAGAAGACTC 
CGTAGATGGGAGTCCGATAT 


(CTDn 


FAM 


57 


122-176 


0.1 


3 


KJ446961 


Usl4 


F: 
R: 


GACCATTAGCCAGGTTAGAG 
CACAGGGTCTCATGATTCTG 


(GCTT) 6 


FAM 


59 


189-205 


0.1 


1 


KJ446962 



Note: T a = annealing temperature. 
a See Appendix 1 for specimens used. 

b Multiplex indicates loci that were mixed in the same capillary electrophoresis run. 
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Table 2 . Characteristics of nine polymorphic microsatellite loci developed 
for Usnea subfloridana* and screened in 174 individuals. 



Locus 


Total ( n 


— y-r) 


EPV ( 


n = 56) 


EVO ( 


n = 60) 


EVY (n 


- SQ\ 

— jo; 


A 


H e 


A 


H e 


A 


H e 


A 


H e 


UsOl 


15 


0.81 


9 


0.81 


14 


0.84 


9 


0.74 


Us02 


6 


0.78 


6 


0.75 


5 


0.79 


5 


0.76 


Us03 


14 


0.78 


9 


0.73 


10 


0.77 


12 


0.83 


Us04 


3 


0.43 


3 


0.48 


3 


0.41 


3 


0.39 


Us05 


9 


0.63 


7 


0.69 


7 


0.63 


6 


0.55 


Us06 


7 


0.71 


6 


0.72 


6 


0.72 


6 


0.71 


Us07 


6 


0.66 


4 


0.65 


6 


0.68 


4 


0.67 


Us08 


7 


0.72 


7 


0.77 


6 


0.70 


6 


0.68 


Us09 


12 


0.42 


6 


0.43 


9 


0.43 


5 


0.41 


Average 


8.78 


0.66 


6.33 


0.67 


7.33 


0.66 


6.22 


0.64 



Note: A = number of alleles; H e = Nei's unbiased gene diversity; n = total 
number of samples analyzed. 

a See Appendix 1 for populations used in the study. 



used to study the genetic differentiation and diversity in U. sub- 
floridana, U. florida, and related species and will allow us to 
investigate effects of forest management and environmental 
pollution on genetic population structure in epiphytic lichens. 
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Appendix 1 . Location information, number of individuals sampled, and voucher specimens of Usnea species used in this study. 



Species 


Locality 


Population 


Latitude 


Longitude 


n 


Herbarium a 


Date 


U. florida * 


United Kingdom: Devon, 


UK 


50°33'20"N 


04°02'58"W 


1 


TU-FLO 02 


07.10.2006 




Dartmoor, Lydford 














U. florida * 


United Kingdom: Devon, Morwellham 


UK 


50°33'20"N 


04°02'58"W 


1 


TU-FLO 03 


06.10.2006 




Quay, river Tamar 














U. glabrescens 


Estonia: Community of Polva Vald 


EPV_16a 


58°06'13.6"N 


27°04'29.1"E 


1 


WSL-21566 


29.06.2007 


U. subfloridana* 


Norway: Hordaland, Voss, Tvinde 


NOR 


60°43'43"N 


06°29'13"E 


1 


TU-SUB 01 


12.07.2006 


U. subfloridana* 


Lithuania: Birzai district, Birzai forest, 


LIT 


56°14'N 


24°58'E 


1 


TU-SUB 05 


13.05.2006 




Buginiai forestry 














U. subfloridana* 


Finland: Northern Savonia, Suonenjoki, 


FIN 


62°44 , 39 ,, N 


27°01 , 28 ,, E 


1 


TU-SUB 10 


10.07.2007 




Mustolanmaki 














U. subfloridana 


Estonia: Community of Polva Vald 


EPV 


58°06'13.6"N 


27°04 , 29.1 ,, E 


56 


WSL-21521-21583 


29.06.2007 


U. subfloridana 


Estonia: Community of 


EVO 


58°07'23.7"N 


26°59 , 22.0 ,, E 


60 


WSL-21773-21834 


03.07.2007 




Vastse-Kuuste Vald 














U. subfloridana 


Estonia: Community of 


EVY 


58°08'29.3"N 


27°02'57.9"E 


58 


WSL-2 1899-2 1961 


04.07.2007 




Vastse-Kuuste Vald 














U. wasmuthii 


Estonia: Community of Polva Vald 


EPV_7a 


58°06'13.6"n 


27°04'29.1"E 


1 


WSL-21539 


29.06.2007 



Note: n = number of samples analyzed; Herbarium = Herbarium accession number, Date = collection date. 

a TU = University of Tartu, WSL = Swiss Federal Research Institute WSL. Specimens stored at WSL are kept frozen at -20°C. All specimens were 
collected by T.T. and L.M. 

* Specimens used for shotgun sequencing. 
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